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Abstract. We investigate the modification of discrete scale invariance in the bound state spectrum by long- 
range interactions. This problem is relevant for effective field theory descriptions of nuclear cluster states 
and manifestations of the Efimov effect in nuclei. As a model system, we choose a one dimensional inverse 
square potential supplemented with a long-range Coulomb interaction. We study the renormalization and 
bound-state spectrum of the system as a function of the Coulomb interaction strength. Our results indicate, 
that the counterterm required to renormalize the inverse square potential alone is sufficient to renormalize 
the full problem. However, the breaking of the discrete scale invariance through the Coulomb interaction 
leads to a modified bound state spectrum. The shallow bound states are strongly influenced by the Coulomb 
interaction while the deep bound states are dominated by the inverse square potential. 
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[ 1 Introduction 

' The application of effective field theory (EFT) methods 
to nuclear physics is by now well established [1] . If there 
is a separation of scales in a physical system, effective 
field theory allows for controlled calculations of low-energy 

■ observables with well-defined error estimates. In nuclear 

■ physics, mostly effective field theories with nucleons and 
pions (and possibly Deltas) as degrees of freedom are used. 

■ However, for a certain class of systems, it is possible to use 
effective field theories with even more effective degrees of 
freedom [2]. This is the realm of halo nuclei and nuclear 
cluster states. 

A halo nucleus is one that consists of a tightly bound 
core surrounded by one or more loosely bound valence nu- 
cleons. The valence nucleons are characterized by a very 
low separation energy compared to those in the core. As a 
consequence, the radius of the halo nucleus is large com- 
pared to the radius of the core. The separation of scales in 
halo nuclei leads to universal properties that are insensi- 
tive to the structure of the core (see, e.g., Ref. [3] and refer- 
ences therein). The most carefully studied Borromean halo 
nuclei are ^He and ^^Li, which have two weakly bound va- 
lence neutrons [HIS] . In the case of ^He, the core is an alpha 
particle. An EFT framework to describe halo systems was 
introduced in [6l[7] where the neutron-alpha (na) system 
was studied in an EFT with nucleon and alpha degrees 
of freedom. Further extensions to the proton-alpha and 
alpha-alpha systems were considered in Refs. [8l[9]. Simi- 



lar concepts can be applied to nuclear cluster states. The 
best-known example is the structure of ^^C. This system 
has an excited 0+ state, the so-called Hoyle state which 
shows a clear clustering into three a particles. This ob- 
servation suggests that this state can be described by an 
EFT of a particles inteacting via short-range contact in- 
teractions. An important question is whether there is an 
universal binding mechanism for these systems. A prime 
candidate for such a mechanism is the Efimov effect [10] . 

In an EFT framework, the Efimov effect can be related 
to a renormalization group (RG) limit cycle [11] . Most ap- 
plications of the RG involve a flow towards a fixed point, 
where the system is scale invariant. However, as pointed 
out by Wilson 12J, one can also have closed curves un- 
der the RG flow in the space of coupling constants. The 
RG flow completes a cycle around the curve every time 
the cutoff is changed by a multiplicative factor Aq. This 
number Aq is the preferred scaling factor. A necessary con- 
dition for a limit cycle is invariance under discrete scale 
transformations: x — *■ XqX, where n is an integer. This 
discrete scaling symmetry is reflected in log-periodic be- 
havior of physical observables. The Efimov effect can be 
understood as the manifestation of a limit cycle in the 
bound state spectrum of the three-body problem with 
large scattering length a. This limit cycle property is man- 
ifest in the EFT treatment of Refs. [13] , where an explicit 
log-periodic three-body counterterm is introduced. In the 
limit a — > ±oo, there is an accumulation of 3-body bound 
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states near threshold with binding energies differing by 
multiplicative factors of Ag ~ 515.03 JJj. Recently, the 
first convincing experimental evidence for this effect was 
obtained by measuring its effect on three-body loss rates 
in a gas of cold Cs atoms [15] . 

The Efimov effect could also be responsible for the 
binding of certain halo nuclei and cluster states 16J. In 
particular for the latter, however, the short-range strong 
interaction is usually accompanied by a long-range Cou- 
lomb interaction. The effect of such long-range Coulomb 
interactions on the physics of limit cycles and discrete 
scale invariance is therefore an important issue. In order 
to get some insight into this question, we start with a 
simpler problem that has also a limit cycle behavior: the 
one-dimensional Schrodinger equation with an attractive 
inverse square potential. If the attraction is larger than a 
certain critical value, the system also shows a limit cycle. 
This limit cycle becomes evident in a bound state spec- 
trum with discrete scale invariance similar to the Efimov 
effect. Indeed, the inverse square potential is intimitely 
connected to the three-body system with large scattering 
length, which reduces to a one-dimensional Schrodinger 
equation with an inverse square potential in the hyper- 
radius for large momenta [Mp]. This makes the inverse 
square potential an ideal model system to study the physics 
of limit cycles and discrete scale invariance pT l lTS l lTQ j lS m 

HUES]. 

In this paper, we study the effect of a long-range Cou- 
lomb interaction on discrete scale invariance in the bound 
state spectrum for an inverse square potential with a Cou- 
lomb potential. (For an earlier study of the interplay be- 
tween Coulomb and strong interactions in exotic atoms, 
see Ref. [2^.) In the next section, we briefly review the 
renormalization of the inverse square potential in the ap- 
proach of Ref. [21]. In Sec. [S] we introduce the long-range 
Coulomb potential. The renormalization and our results 
for the bound state spectrum are discussed in Sec. [H and 
in Sec.[5|a perturbative treatment of the Coulomb interac- 
tion is given. Finally, we present our conclusions in Sec.[6j 
Our treatment of the Coulomb divergence is described in 
Appendix [XJ 



2 Inverse Square Potential 

In order to set up our problem, we briefly review the renor- 
malization of the l/r"^ potential in momentum space with- 
out the long-range Coulomb interaction ^21j . We consider 
the attractive inverse square potential 



Vs{r) = 



with r = Irl 



_ 1 2 



(1) 



and 1/ > a positive real parameter. This potential has the 
same scaling behavior as the kinetic energy operator and, 
consequently, is scale invariant at the classical level. In the 
following, we set the particle mass and Planck's constant 
m = ft = 1 for convenience. For values of c > — |, the po- 
tential is well-behaved and the corresponding Schrodinger 
equation has a unique solution, see Ref. [241. However, 



we are interested in the case c < — ^ which corresponds 
to real values of in ([1]). In this case, the l/r^ poten- 
tial is singular and the usual boundary conditions for the 
Schrodinger equation do not lead to a unique solution. We 
can calculate the Fourier transform of the potential using 
dimensional regularization (see Ref. ^2lj for details). This 
leads to the expression 



Vs{q) 



q 



(2) 



for the momentum space representation of the 1 /r'^ poten- 
tial ([TJ. Since the potential is local, its Fourier transform 
depends only on the momentum transfer q. 

The Lippmann-Schwinger (LS) equation for two parti- 
cles interacting via Vs{q) from Eq. ([2]) in their center-of- 
mass frame takes the form 



T£(p,p') = ^s(|p-P'| 



rf'g ^5(|p-q|) 
(27r)3 E -q^ +ie 



Te{ci,p'), 
(3) 

where E is the total energy and p (p') are the relative mo- 
menta of the incoming (outgoing) particles, respectively. A 
pictorial representation of this equation is given in Fig. [1] 
We are only interested in the S-wave contribution. In higher 





Fig. 1. Lippmann-Schwinger equation for a two-body potential 
V. 



partial waves the singular behavior of the potential is 
screened by the angular momentum barrier, but for suf- 
ficiently strong attraction it will become visible as well 
(see, e.g., Ref. [22] )• Projecting onto S- waves by integrat- 
ing the equation over the relative angle between p and p': 
J d cos 9pp' , we obtain the integral equation 



tE{p,p') = vs{p,p') + 



27r2 



dq q^ 



E-q^ 



X vs{p, q) tE{q,p') : 



where 



vsiP, q) = 271-^0 



r ejp-q) ^ e{q-p) 
VP q 



(4) 



(5) 



The physical observables are the bound state spec- 
trum and the scattering phase shifts 6{k). The phase shifts 
are determined by the solution to Eq. ^ evaluated at 
the on-shell point E = fc^, k = p' — p via ti^2{k,k) = 
—Att/ (fc cot 5{k) — ik) . Since p' appears only as a parame- 
ter in Eq. ([?]), we can set p' = p to simplify the equation. 
The binding energies are given by those values of -E < 
for which the homogeneous version of Eq. ^ has a solu- 
tion. For the bound state equation the dependence of the 
solution <j)E{p) on p' disappears altogether. 
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It is well-known that Eq. (j4|) does not have a unique 
solution since the Xjr^ potential for real v is singular [24 . 
The most general solution of the bound state equation for 
i? = can be written as 



00 (p) ^Mv 



-1/2 (pi^e^V 



■V 



-lip 



(6) 



where the relative phase is a free parameter. The value 
of if is not determined by the l/r^ potential and has to 
be taken from elsewhere. It is exactly this phase which 
is fixed by self-adjoint extensions of the potential [18]. 

In the framework of an effective theory, this is conve- 
niently done using renormalization theory. We regularize 
the LS equation by applying a momentum cutoff A and 
include a momentum- independent counterterm 5Vs{A) in 
the potential. The precise form of the cutoff, for example 
Gaussian cutoff or sharp cutoff, is not important, but we 
use a sharp cutoff for simplicity. Making the replacement 



Vs{q) 



Vs{q) + 5Vs{A) = 



27r2c Hs{A) 



q 



A 



(7) 



in Eq. ^ the LS equation (jH) for bound state solutions 
with E = —Eb < becomes 



<{P) = - 



dqq^ 



1 



2^2/0 EB+q" 



vs{p, q) 



Hs{A) 



A 



(8) 



The functional dependence of Hs{A) can be determined 
analytically from invariance of low-energy observables un- 
der renormalization group transformations. We demand 
that the relative phase of the zero-energy bound solution 
of Eq. ([HI) remains unchanged under variations of the cut- 
off A and find 



Hs{A) = 2tt\: 



2 1 - 2i^ tan(i^ ln(yl/yl,)) 



1 + 2vt(m{v\n{A/A^)) ' 



(9) 



where vl* is a free parameter that determines the relative 
phase in ([6]): = —vlnA^,. In order to fix ip, we can either 
specify both the cutoff A and the dimensionless coupling 
H or, using Eq. ([9|), one dimensionful parameter: A^,. This 
parameter yl* is generated by the iteration of quantum 
corrections in solving the integral equation ([H]). This is 
similar to the phenomenon of dimensional transmutation 
in QCD pg. 

Note that Hs{A) remains unchanged when the argu- 
ment is multiplied by Ag, where n is an integer number 
and Ao — e^l^ is the discrete scaling factor. This discrete 
scaling symmetry is a consequence of the limit cycle and 
reflects itself in physical observables [2jl2T| . In the bound 
state spectrum, for example, the ratio of consecutive bind- 
ing energies is a constant, En/En+i = . Therefore, 
we can study modifications of the limit cycle through mod- 
ifications to the discrete scaling symmetry of the bound 
state spectrum. 

Moreover, the discrete symmetry implies the existence 
of a set of cutoffs 



An{A^) — yl* exp(n7r/z^) 



(10) 



with i/5(yl„) = 0. We can therefore obtain a renormal- 
ized version of Eq. ([8]) that does not explicitly contain 
the counterterm by using the discrete set of cutoffs from 
Eq. (fTO|l . The same trick can be used for the three-body 
problem with large scattering length 26J. 



3 Inclusion of the Long-Range Potential 

We now include an additional attractive Coulomb poten- 
tial of the form 



yc{q) = - 



q^ -\- 7^ 



(11) 



where 7 is a photon mass that will be taken to zero in 
the end. The speed of light has been set to unity for con- 
venience. In the following, we vary the strength a of the 
potential. Projecting onto S-waves as discussed in the pre- 
vious section we find 



vc{p,q) = -—In 

pq 



{P + q? 



-7 



{p - qY + 7^ 



(12) 



The integral equation for boimd state solutions, Eq. ([8|), 
then becomes 



4>Eb (p) 



1 

'2^ 



dqq^ 



EB + q^ 
Hs{A) 



vs{p, q) 



A 



+ «c(p, q) 



4>EAq)- (13) 



From Eq. p^ . it is clear that the kernel of the inte- 
gral equation (I13|) diverges for q = p m. the limit 7 — > 0: 
this is the well-known Coulomb singularity. Its origin can 
be traced back to the integral equation for the scattering 
amplitude with a long-range Coulomb interaction, which 
diverges at forward angles. When projected into S-waves, 
this singularity appears in the diagonal terms of the po- 
tential in momentum space. For the binding energies this 
singularity should not be a problem, as long as the di- 
agonal terms are handled properly. In Appendix [Al we 
describe how to treat these terms, based on the the idea 
outlined in Ref. [23. 



4 Renormalization and Results 

First we address the renormalization of the full problem 
including Coulomb. It is not clear a priori whether the 
introduction of the Coulomb potential will require an ad- 
ditional counterterm. 

In order to answer this question, we calculate the bound 
state spectrum at a given cutoff and choose to fix the bind- 
ing energy of one bound state as the cutoff A is varied. The 
fixed energy was chosen as one of the pure l/r^ states, 
^'fixcd ~ 0.3080 We then calculate the bound state spec- 
trum for other values of the cutoff A and adjust the coun- 
terterm Hs{A) numerically to keep this binding energy 



^ With our choice of units all energies and momenta are di- 
mensionless. 
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Fig. 2. Bound state spectrum with one state fixed at Eb ~ 0.308 (left panel) and the counterterm required for renormalization 
(right panel) for a = 0, 0.1, 0.3, 0.5 and with = yl* = 1, as a function of A. In the right panel, the curve a = coincides with 
H{A) given by Eq. Q, except for /I < 1. 



fixed. The results of this analysis are shown in Fig. [51 The 
left panel shows the bound states for a — 0, 0.1, 0.3, 0.5 as 
a function of A. In order to keep this and the remaining 
figures as legible as possible, only the three deepest states 
are shown. The case a = corresponds to the pure 1/r^ 
potential. One observes that the cutoff dependence of all 
states is removed once the counterterm is adjusted to fix 
one of the states. Note, however, that the deepest state 
shows a cutoff dependence near the cutoff where it first 
appears with infinite binding energy. This behavior has 
nothing to do with the long-range interaction and is due 
to the way the system is renormalized, keeping low-energy 
physics unchanged. A similar behavior is also observed for 
the pure 1/r^ potential and the three-body system with 
large scattering length ^SjHHEI]- 

Moreover, it is evident that the long-range interac- 
tion destroys the discrete scale invariance in the spectrum. 
Only for a = 0, the ratio of consecutive binding energies 
is a constant. In the right panel of Fig. [21 we show the 
numerically determined values of the counterterm Hs{A) 
as a function of A for a — 0,0.1,0.3,0.5. The figure sug- 
gests that the long range potential shifts the argument of 
the counterterm ^^^(Inyl) Hs{\nA + f{a)) where f{a) 
is a monotonic function of the coupling a. This implies 
that the long-range Coulomb potential merely renormal- 
izes the value of A^. As a consequence, the counterterm 
Hs{A) from Eq. ([9]) should be sufficient to renormalize 
Eq. (fT3|) including the long-range Coulomb potential. 

In order to test this assumption, we calculate the bound 
state spectrum using the counterterm for the pure l/r^ 
potential from Eq. ([9]). The result is shown in Fig.[3l For 
simplicity, we take v = A^ = 1 and show different values 
of a. For each energy level, the binding energies increase 
as the Coulomb strength is increased. It is evident that the 
Coulomb potential infiuences the spectrum, but does not 
destroy the cutoff independence of the binding energies. 
Clearly, no additional counterterm is required. 

Note also that the shift in energy due to the Coulomb 
interaction depends on the exitation level — the deepest 
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Fig. 3. Binding energies for f 
cutoff yl for a = 0, 0.1, 0.3, 0.5. 
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bound states have the smallest shifts relative to the pure 
l/r^ case, while the shallower ones have larger shifts and 
resemble more the Coulomb spectrum. This behavior can 
be understood by comparing the relative strengths of the 
l/r^ and l/r potentials in coordinate space. There will 



always be a special distance r — \c\/a = (1/4 



where the two potentials have equal strength: Vs{f) 
Vc{r)- The qualitative pattern of the energy levels can be 
understood by comparing the binding energy of a given 
state, Eb, with the potential energy at f, 



E, 



pot 



\Vsif) + Vcif)\ 



2a' 



1/4- 



(14) 



The bound states with Eb S> -Epot are more sensitive to 
shorter distances where Vs dominates over Vc- Therefore, 
the spectrum resembles the 1 /r^ spectrum and is approxi- 
mately scale invariant. For Eb ^ Epot, on the other hand, 
the states are more sensitive to larger distances and the 
spectrum is similar to the Coulomb spectrum. 
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An alternative way of understanding this behavior is 
by looking at Eq. (jl3|) . We take one of the bound states 
Eb = Bn+i with an eigensolution (f)n+i{p). The next 
deeper bound state, denoted by i3„, defines a variable i„ 
(> 1) via t"^ = Bn/Bn+i- Bn has an eigensolution 
that satisfies 



i{p) = -77-2 



dq 



B„ 
Hs{A/A,) 
A 



9{p-q) , eiq-p) 



na 

— In 

pq 



P 

(p+q)^ 



ip-qf 



Mq), (15) 



where we explicitly indicated the dependence of Hs on 
A and yl*. Rescaling the external momentum p and the 
integration variable q by t„ yields 



dqq'^ 



I 



,/o s„A2+r 



2tt^c 
1 na 



0{p-q) , 0{q^p) 



HsitnA/A, 



A 



■In 



tn pq [{p-q) 



{p+qf 



(t)n{tnq)- 



In addition, we set A tnA assuming cutoff indepen- 
dence, which is verified numerically. Writing tn = exp(7r/iy+ 
6n) and using the log-periodicity of Hs, we conclude from 
Eq. (fT6|) that -B„+i is also a bound state of a Hamiltonian 
that has a Coulomb potential weakened by l/i„ and pa- 
rameter A^ multiplied by e~^^ . By induction, it follows 
that a deeper state Bn-k, with eigensolution (pn-k{p), 
is given by akBn+i = {n^^Qtn-j)Bn+i. Furthermore, 
Bn+i is an eigenvalue, with eigensolution 0„_fc(crfcp), of 
the Hamiltonian Hq + Vs + Vc jok with a parameter A^ 
multiplied by e~^'' = exp(— ^^~q (5„^j). Therefore, the 
solution (pn-kip) with eigenvalue Bn-k is only weakly sen- 
sitive to the Coulomb potential. The same rationale ap- 
plies to the solution (j)n-k-i{p) for the next deeper state 
Bn-k-i- One therefore expects that Sn-k-i tends to zero 
and, consequently, tn-k-i to the discrete scaling factor 

In order to support these conclusions, we plot in Fig.|4] 
the binding energies as functions of the "excitation num- 
ber" n, the latter being defined relative to the deepest 
bound sta,te at yl = 2000 For = 1 and a = 0.1,0.3,0.5, 
one has Epot = 1.6 x 10-^ 1.44 x 10-\4.0 x 10-\ re- 
spectively. The exact scale invariance is broken for any 
finite value of a. The figure illustrates how the limit of 
exact discrete scale invariance is approached as the states 
become deeper. It confirms that for Eb ^ ^pot the behav- 
ior is closer to the geometric 1/r^ spectrum with In Eb = 
const. x2tt/i^. This spectrurn is indicated by the solid 
straight line. For large n, where Eb <^ -Epot, the spectrum 
approaches the Coulomb spectrum. This is illustrated by 



^ Note that this assignment is not unique since new deep 
states appear as the cutoff is increased. In the EFT frame- 
work, this is not a problem since all states outside the range 
of validity of the EFT can be ignored. 
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Fig. 4. Binding energies for 1/ = At, — 1 and a = 0, 0.1, 0.3, 0.5 
as function of the "exitation number" n (see text for details). 



(16) the dotted lines which represent the Coulomb energies 
Q;^/(4n'^), with the shifted excitation number n = n — 2. 
This particular choice of h is natural since the n — S level 
is the closest to the n = 1 level of the pure-Coulomb spec- 
trum. In accordance with our expectations, the Coulomb 
pattern is already observed at moderate n for a = 0.5 and 
0.3, while it is achieved at larger n for a = 0.1. 



Next we study the dependence of the bound state spec- 
trum on v and A^ . In Fig. [5l we show the spectra for v = 1 
and A^, = 2 (left panel) and for yl* = 1 and f = 2 (right 
panel). We first consider the left panel. Due to the dis- 
crete scale invariance, a change of A^ in the pure 1/r^ 
case modifies the values of the energies but preserves the 
geometric character of the spectrum, i.e. the ratio of subse- 
quent states is determined by the preferred scaling factor 
squared (e'^/'^)^. This behavior is manifest as a vertical 
displacement of the spectrum a = relative to the one 
shown in Fig. [3l A similar shift of the energies is also 
expected for a ^ 0. However, the Coulomb interaction 
breaks the discrete scale invariance and modifies the ra- 
tios between consecutive energy levels. The increase of vl* 
provides more binding to the system, while the splitting 
among different values of a is reduced. This observation 
is in agreement with the finding in the previous section: 
deeper bound states become less sensitive to the Coulomb 
part of the interaction, which is responsible for the dif- 
ferent splittings. Next we consider the case yl* = 1 and 
1^ — 2 shown in the right panel of Fig. [5] For a given value 
of a, the 1/r'^ part of the interaction becomes relatively 
stronger and the overall binding is increased. The influ- 
ence of the long-range Coulomb part is decreased. For the 
pure l/r^ case a = 0, the discrete scaling factor is reduced 
to e'^/^ and the states move closer together. For the case 
a ^ 0, the spectrum has distinct features of the l/r^ prob- 
lem — very little spread in energy among the different a 
considered, and a binding ratio very close to (e'^/^)^. 
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Fig. 5. Bound state spectrum for /I* = 2 and v = \ (left panel) and for A^, = \ and v = 2 (right panel). 
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5 Perturbative Treatment of the Coulomb 
Part 

In the previous section we saw that the deepest bound 
states have small Coulomb corrections. In this section we 
calculate these energy shifts relative to the unperturbed 
solution in perturbation theory in order to test our 
hypothesis. 

It is convenient to start from the LS equation for the 
scattering state instead of the transition amplitude ([3]). 
The former reads 



) ^ |p) + Gf\E) [Vs + AVb) , (17) 



homogeneous version of Eq. p9|) to eliminate the state 
|(/?p~''^) from the above equation. This is achieved via mul- 
tiplication from the left by (xp' | Vs and yields 



;(+)/r(0), 



(21) 



We can again use Eq. to express e'^^ = ~E^^^ in 
terms of the transition amplitude Te{p,p') defined by 
Eq. ([3]). After integration over the angles, we obtain 



= ~A/B 



(22) 



where Hl'f)) is the initial (final) state, G^o^'{E) 

is the free two-particle propagator, and A is a parameter 
expansion that will be set to 1 in the end. Expanding the 
scattering state and the binding energy in powers of a 
small parameter A, 



,(-) 



E^E^°'> +XE^^^ + ■ 



(18) 



one gets to leading order in A the LS equation for the pure 
l/r^ interaction Vg: 



lx^+^> = |P> 



(19) 



Multiplication of this equation by (p'|Vs gives the LS 

equation for Te{p,p') = (p'|Vs|x^+^) (cf. Eq. Col- 
lecting the terms linear in A, we obtain 



Gi+\E('y) Vc - E^'^ G^+Hi?(°))2 Vs 
+G(+)(i?("))Vs|^(+V 



(20) 



This is an integral equation for the state Iv'p^'') with a 
driving term proportional to the unperturbed state |xp^'). 
Since we are dealing with bound states, one can use the 



where 

A = 



A 

dp 

.^0 



dq^Mv. r"C'(p, q) 

B 



B = 27r' 



"P — i7r\ 

{e'^b +P^)^ 



g'0(O)(g) 
E 



A 



(23) 



and (/)(°)(p) is the solution of Eq. 

A comparison of the exact binding energies with the 
perturbative results is given in Fig.[6l We show the pertur- 
bative Coulomb binding energies for a = 0.1 (left panel) 
and 0.5 (right panel), compared with exact energies and 
the energies for the pure 1/r^ case. On one hand one 
clearly sees that the perturbative treatment works quite 
well for the two deepest bound states, where the effect of 
the Coulomb interaction is expected to be small. This is 
true even for a relatively strong Coulomb potential with 
a = 0.5. The shallowest state in Fig. \6\ on the other 
hand, cannot be described by perturbation theory in the 
Coulomb potential. In this case, we no longer have Eb ^ 
^pot (cf. Eq. (fT4|) ) and the Coulomb effects are large. In- 
deed, the perturbative treatment of the Coulomb interac- 
tion gives 4.4 X 10~^ (1.9 x 10~^) compared to the exact 
value 6.3 x 10"^ (5.2 x lO'^), for a = 0.1 {a = 0.5). These 
results clearly support our hypothesis from the previous 
section. 
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Fig. 6. Bound state spectrum witli a perturbative Coulomb treatment for a = 0.1 (left panel) and a 
compared against the exact and a = cases. The parameters At and v were set to 1. 



100 

0.5 (right panel), 



We should point out that a perturbative treatment of 
the l/r^ potential relative to the Coulomb potential to 
calculate the shallower states is not possible. This is due 
to the singular nature of the 1 /r^ potential for real values 
of V in Eq. ([T]) which we consider here. We have verified 
explicitly that if the potential is not singular (cor- 
responding to imaginary v) the perturbative treatment 
works quite well. The latter case, however, corresponds to 
a situation where the Schrodinger equation has a unique 
solution and limit cycles are absent. 

6 Summary and Conclusions 

In this work, we have investigated the modification of limit 
cycles and discrete scale invariance by the presence of a 
long-range interaction. As a specific example, we have con- 
sidered the quantum mechanical inverse square potential 
supplemented by an attractive long-range Coulomb inter- 
action. We have focused on the bound state properties of 
this model system. 

Our study of the cutoff dependence of the binding en- 
ergies shows that no additional counterterm is required 
for renormalization when the Coulomb potential is added. 
The counterterm that renormalizes the inverse square po- 
tential alone is sufhcient to renormalize the full problem. 

In the presence of the Coulomb potential, the countert- 
erm can no longer be obtained analytically. We have cal- 
culated the counterterm numerically by fixing one of the 
bound state energies. All other binding energies are then 
independent of the ultraviolet cutoff A. This procedure has 
been carried out for several values of Coulomb strength 
parameter a. The counterterm was confirmed to be a log- 
periodic function with discontinuities. Its yl-dependence 
is the same as for the pure inverse square potential but 
shifted along the A-asds. Such a translation corresponds 
to a finite renormalization of the counterterm parameter 
/I*. 

The discrete scale invariance of the inverse square po- 
tential is broken by the Coulomb potential. We have inves- 
tigated the deviations from discrete scaling symmetry for 



different strengths of the Coulomb potential. For highly 
excited states, the long-distance Coulomb tail dominates 
the dynamics and the levels tend towards a Coulomb spec- 
trum. The deepest bound states, however, are mostly sen- 
sitive to the short-range potential and show an ap- 
proximate scaling symmetry. The spectra obtained for var- 
ious Coulomb strengths and values for A^, and v were stud- 
ied in detail. Due to the breaking of scaling symmetry, the 
ratio between consecutive bound state energies (the dis- 
crete scaling factor) is no longer a constant. The splittings 
of the energy levels depend on the magnitude of the bind- 
ing and the Coulomb strength. 

To verify our conclusions, we have studied the behav- 
ior of the deep bound states in perturbation theory. We 
have derived an expression for the energy shift relative to 
the spectrum, treating the Coulomb interaction in 

first-order perturbation theory. We have shown that the 
perturbative expression is valid for bound states that sat- 
isfy Eb ^ -Epot, where -Epot is the potential energy at the 
distance f where both interactions have equal strength. 

After the general features of the breaking of discrete 
scale invariance are understood for this example, we are in 
the position to study more realistic systems. Our results 
could be useful for the study of nuclear cluster states in 
the halo EFT [51[71[51[g] . In Ref. e.g., a power counting 
scenario for the aa system was formulated. According to 
this scenario the *Be system would exhibit conformal in- 
variance at leading order and ^^C would display an exact 
Efimov spectrum. These exact features are broken by the 
Coulomb interaction but some remnants of this behav- 
ior are manifest in the experimental spectra, such as the 
shallowness of the ^Be 0+ resonance. The ^^C Hoyle state 
would then be a remnant of an Efimov state that appears 
in the limit of large scattering length. An application of 
the power counting scenario to the triple-a system re- 
mains to be carried out. Our calculation provides a first 
step towards the understanding of the breaking of discrete 
scale invariance in these systems. Additional expansions 
such as a strong coupling expansion for the Coulomb in- 
teraction ^ might be useful and deserve further study. 
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A Treatment of the Coulomb Divergence 

In this Appendix, we describe our method to treat the 
Coulomb singularity in Eq. p^ . 

The integral equation for bound states ([13]) is conve- 
niently rewritten afl 



vsip, q) 



+vc{p,q) (pE{q) ~ F{p)(j)E{p) , 

(24) 

where the function F{p) was introduced on both sides of 
the equation to cancel the Coulomb divergence in the di- 
agonal terms. Following Ref. [27], one finds that a suitable 
choice for F{p) is given by 



1 

2^ 



"-11 — vc{p,q) . 



/o {E-q^Y 
The integral can be evaluated analytically: 

a{E- p2) 



(25) 



Fip) - 



(i;-p2+^2)2+4p2^2 

f {E-p'^-^i'^'^ 



arctan 



A 



arctan 
1 



A + p 



arctan 



P 



A'Kp{E-A^) 



1 J \ 1 

{E-p'W){A'-p'W) 



+Ap^-i'^ 



In 



[p+Af + "f 



21 



(p-yl)2 + 72 
a{E-p^) 



(26) 



[E^p-^ 



v2^2 



4j)2^2 



_ p2 _ ^2 



+ 7 



(27) 



Using Eq. ([25]) on the r.h.s. and either Eqs. ([26]) or ([27]) 
on the l.h.s. of Eq. ([24]l is enough to eliminate the diver- 
gence problem. As a numerical check, we set and 
obtained very stable and accurate values for the Coulomb 
spectrum. For the values of the Coulomb strength a con- 
sidered in this work, we also observed cutoff independence 
except at lower values {A < 1), where results using (|26p 
or ([77]) start to deviate by a few percent. 



^ We do not explicitly display the counterterm Hs{A) since 
it is irrelevant for the Coulomb divergence problem. 
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